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MAIN  BATTLE  TANK  FLEXIBLE  GUN  TUBE  DISTURBANCE  MODEL 
THREE  SEGMENT  MODEL 


Dr.  Henry  J.  Sneck1 

'U.S.  Army  Benet  Labs,  TACOM-ARDEC,  Building  #40,  Watervliet,  NY  12189-4050 


A rational  approach  to  disturbance  rejection  is  proposed  and 
applied  to  a simple  three-degree-of-freedom  flexible  gun  tube 
model  using  feedforward  and  feedback  compensation.  The  first 
two  natural  frequencies  of  the  pin-free  and  cantilever  tube  are 
matched  by  adjusting  the  dimensions  of  the  rigid  segments  and 
the  stiffness  of  the  torsional  springs  that  join  them.  It  was  found 
that,  contrary  to  the  previously  analyzed  two  degree-of-freedom 
segment  model,  the  muzzle-end  segment  could  be  stabilized  by 
the  proper  choice  of  transfer  functions  and  elevation  driveline 
response.  The  analysis  serves  to  establish  the  requirements  for 
the  transfer  functions  and  stabilizing  actuator  systems. 


INTRODUCTION 

Modern  tank  cannon  are  long,  relatively,  thin,  beam-like  hollow  cylinders.  Their 
accuracy  is,  in  part,  determined  by  their  flexibility,  especially  under  dynamic  loading.  Very 
small  deflections  and  rotations  of  the  muzzle  end  can  have  a significant  influence  on  the 
accuracy  of  the  shot  at  long  ranges.  Muzzle  motions  induced  by  firing  are  inevitable,  and 
difficult  to  control  because  of  the  time  scale  of  the  firing  is  of  the  order  of  milliseconds. 

Another  source  of  muzzle  motion  is  the  ground-induced  motion  of  the  vehicle.  These 
motions,  transmitted  through  the  trunnions  and  gun  actuators,  can  be  quite  large  and  have 
frequencies  comparable  to  the  natural  frequencies  of  the  tube.  The  time  scales  of  these 
disturbances  depend  on  the  tank  speed  and  on  the  nature  of  the  terrain.  They  are  typically  of  the 
order  of  seconds  or  longer.  Sensing  and  actuation  to  control  the  influence  of  vehicle  motion  on 
the  muzzle  response  might  be  possible,  given  these  relatively  long  time  scales.  This  raises  two 
questions.  First,  is  it  possible  to  reject  some,  or  all  he  ground  motion  disturbance  from  the 
muzzle  motion?  In  a previous  paper  [1]  it  was  suggested  that  not  all  of  the  disturbance  could  be 
rejected.  Second,  if  the  more  comprehensive  model  used  here  indicates  that  all  of  the 
disturbance  can  be  rejected,  what  is  the  required  control  strategy? 

During  a discussion  with  Dr.  Purdy,  author  of  Reference  [2],  he  suggested  that  the 
fidelity  of  his  two-segment  flexible  model,  documented  in  Reference  [1],  was  inadequate.  He 
recommended  that  the  tube  should  be  divided  into  at  least  three  segments,  with  intervening 
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torsional  springs  and  dampers.  The  author  is  indebted  to  Dr.  Purdy  for  this  suggestion  since  this 
paper  is  the  result  of  his  recommendation. 

EQUATIONS  OF  MOTION 

Figure  1 shows  the  generic  model  of  the  tube  and  the  various  quantities  that  determine  its 
dynamic  behavior. 


A free  body  analysis  of  this  model  yields  the  classical  dynamic  equation, 


M H+H  ^}+[^]  fib  i1)  (!) 

The  elements  of  the  mass,  damping,  and  stiffness  3x3  matrices  are  shown  Appendix  A. 

Purdy  [3]  has  shown  that  the  tube  motion  can  be  adequately  modeled  if  the  segmented 
model  matches  the  pinned-free  and  cantilever  frequencies  of  the  mounted  tube.  Matching  is 
accomplished  by  adjusting  the  size  of  the  rigid  segments  and  the  stiffness’  of  their  connecting 
torsional  springs.  The  2x2  matrices  for  the  cantilever  mode  are  given  in  Appendix  B. 

Transformation  of  Equation  (1)  into  the  frequency  domain  will  allow  its  incorporation 
into  the  control  strategy.  Taking  the  Laplace  transform  of  Equation  (1)  yields: 
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The  elements  of  [a]  and  [I]  are  listed  in  Appendix  C.  This  equation  relates  the  response 
vector  on  the  left  to  the  disturbance  vector  on  the  right.  These  vectors  also  contain  the  actuator 
force,  Fc,  and  the  actuator  displacement,  xc , in  addition  to  disturbances  (s2yt(s),  s0p(s),  0p(s))  and 
the  responses  0i(s),  02(s),  03(s). 
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where  [C]  is  the  transpose  of  the  numerators  of  the  cofactors  of  [a].  The  elements  of  [C]  and  det 
[a]  are  listed  in  Appendix  D. 

The  final  step  in  the  preparation  of  the  dynamic  equations  is  to  perform  the  operation 
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det[a] 


[C]  [/]  = [*] 


(4) 


where  the  elements  of  [B]  can  be  found  in  Appendix  E. 


The  result  of  these  straight  forward,  but  laborious  manipulations,  is  an  equation  for  the 
response  to  the  disturbance  in  terms  of  the  properties  of  the  model  contained  in  [B],  i.e., 
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Of  course  it  was  known  at  the  outset  that  Equation  (1)  could  be  put  into  this  form.  This 
section  merely  provides  the  details  of  how  this  transformation  is  performed,  and  documents  the 
intermediate  steps  and  their  components. 

FEEDBACK  AND  FEEDFORWARD  CONTROL 

The  portion  of  the  response  due  to  the  applied  actuating  force  is 
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where  [Gp]  is  the  "plant"  transfer  function. 


The  portion  of  the  response  due  to  the  disturbance  is 
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where  [Gd]  is  the  disturbance  transfer  function  and  { D } is  the  disturbance  vector. 


Figure  2 is  a block  diagram  of  the  gun  tube  system  with  a gain  Gc,  feedback  H,  feedforward  GCd, 
and  reference  signal,  R.  Because  R is  a scalar  the  feedforward  transfer  function  is  a row  vector, 
i.e. 


[gJ=[ghg12g13] 


(8) 
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Assuming  that  the  tube  rotations  at  the  trunnions  and  the  muzzle  can  be  sensed,  the 
feedback  transfer  function  is  also  a row  vector,  i.e., 

[//]  = [0  G22  0 G24]  (9) 

Now  referring  to  Figure  2 the  response  to  the  disturbance  for  the  controlled  system  is 


10,]  ■ - IgJg,  ] [c,  Jo}= [[/] +[gJ[g„]  [*#„}  (10.) 


or  in  abbreviated  notation 


[d}{D}=[q]{eD} 


(10b) 


The  final  step  in  these  manipulations  is  to  solve  for  the  response  to  the  disturbance,  which  is 

KMd'1  k]M  = ^][>f  \d\{D}  (10c) 


where  [r]  is  the  matrix  of  the  cofactors  of  [qj. 

The  matrices  [r]T  and  [d]  are  given  in  Appendix  F.  It  is  interesting  to  note  that  [d] 
contains  only  the  feedforward  transfer  functions  Gn,  G12,  and  G13,  while  [r]  and  det  [q]  contain 
only  the  feedback  transfer  functions,  G22,  and  G24. 

DISTURBANCE  REJECTION 


Referring  to  Appendix  F the  expanded  version  of  Equation  (10c)  is 
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det [<y ] - 1 + B2]  Gr  G22  + S4I  Gc  G24  (11) 


To  completely  remove  the  effect  of  the  disturbances  on  03(s)d  requires  that 
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and 


r24  d2]  + ru  d4l  0, 
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One  way  to  accomplish  this  is  to  let  G22  = 0 so  that 

r24  = - B4l  G,  G22  = 0 (13) 

and  then  choose 

d4[  = Bn  - Bxl  G(  Gn  =0,  (14a) 

dxl  = B4i  - Bm  G(  Gi2=0,  (14b) 

d4i  = BU-  Bm  Gc  G,  3 = 0 (14c) 

so  that 

det[g]  = 1 + B4I  G,  G24  = rt]  = r22  = rix  (14d) 

and 
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The  effect  of  this  choice  on  the  disturbance  transfer  function  is 
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(16) 


An  alternative  strategy  is  to  let  144  = 1 + B2i  GCG22  = 0,  and  then  choose  d2i  = d22  = d23 
0 so  that  det  [q]  = B41  Gc  G24  = rn  = ^3.  However,  in  the  end  the  resulting  disturbance  transfer 
function  is  the  same  as  created  by  Eq  (16). 
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Since  the  fourth  column  of  [r]T is  eliminated  by  the  matrix  multiplication  and  det[q]  = rn  = r22  = 
r-53,  det[q]  will  be  eliminated  from  the  transformation. 

IMPLEMENTATION 

This  analysis  of  the  three-segment  model  indicates  that  model  muzzle  element  can  be 
stabilized  by  properly  selecting  the  feedforward  and  feedback  transfer  functions.  This  is  contrary 
to  the  finding  for  the  lower  order  two-segment  model  [1],  Although  the  three-segment  model 
only  approximates  the  real  tube,  the  results  of  this  model  are  encouraging  with  respect  to  real 
tubes. 


In  order  to  achieve  muzzle  stabilization  the  breech-end  of  the  gun  must  be  actuated. 
Segments  1 and  2 will  also  rotate.  These  motions  are  determined  by  Equations  (14)  and 
Equation  (16).  All  of  the  elements  of  [d]  can  be  written  in  terms  of  the  elements  of  [B]. 
Because  [r]T  acts  like  an  identity  matrix,  the  product  [r]T[d]  is  quite  simple.  If  the  process  is 
carried  a step  further  the  result  can  be  put  in  terms  of  the  elements  of  [C]  with  startling  results, 
i.e.,  d^2  = d23  = d32=  d33  = 0,  and  only  dn  = d]2  = dj3  = d2i,  and  d3i  are  non-zero.  The  surviving 
elements  are 
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These  transfer  functions  relate  the  disturbances  to  the  responses.  All  of  the  C’s  are  of 
order  s4  with  the  exception  of  Cn,  which  is  of  order  s6.  Since  there  is  no  restraining  torsional 
spring  connecting  the  tube  to  the  mount  in  the  model  s = 0 is  a root  of  det[a].  Removing  this 
rigid  body  factor  from  det[a]  reduces  it  to  order  s5. 

TUBE  MODEL  PARAMETERS 

The  feedforward  transfer  functions  depend  on  the  length  and  mass  properties  of  the 
segments,  the  torsional  stiffness  of  the  joining  springs,  and  the  torsional  damping  coefficients. 
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These  are  chosen  so  that  the  actual  cantilever  and  pin-free  mode  shapes  and  natural  frequencies 
are  matched  as  closely  as  possible  [3].  To  simplify  the  matching  process  it  is  assumed  that  the 
damping  is  negligible.  The  first  estimate  of  the  segment  lengths  can  be  obtained  by  "fitting"  the 
straight-line  segments  to  the  mode  shapes  obtained  from  a finite  element  model  of  the  tube  or 
other  modal  analyses.  This  fitting  is  best  done  by  graphically  overlaying  the  segments  on  plots 
of  the  mode  shapes  to  estimate  the  segment  lengths.  The  calculation  of  the  mass  properties  of 
the  segments  can  then  be  performed  and  these,  along  with  the  modal  frequencies,  inserted  into 
the  characteristic  equations.  The  characteristic  equations  will  then  contain  only  the  torsional 
stiffnesses  as  unknowns.  The  cantilever  and  pin-free  equations  are  both  quadratic  so  that  the 
stiffness  coefficients  can  be  found  directly.  The  degree  of  matching  is  determined  by  how  close 
the  cantilever  and  pin-free  stiffnesses  agree. 

The  characteristic  equations  for  the  cantilever  and  pin-free  segments  are  given  in 
Appendix  G.  Although  the  pin-free  equation  appears  to  be  sixth-order  it  has  a double  root  that  is 
zero.  The  calculations  for  this  trial-and-success  process  are  easily  implemented  on  a spreadsheet. 

The  XM291  tank  gun  was  chosen  for  modeling  because  its  mode  shapes  and  frequencies 
were  available  from  an  existing,  validated  analytical  model.  Matching  the  stiffnesses  proved  to 
be  surprisingly  easy,  requiring  only  modest  adjustments  to  the  first  estimates  of  the  segment 
lengths.  Since  all  their  frequencies  (cantilever:  97.4  Hz,  40.35  Hz;  pin-free:  25.08  Hz,  81.59 
Hz)  were  inserted  into  the  characteristic  equations,  they  are  matched  exactly.  The  torsional 
stiffness  for  the  pinned-free  and  cantilever  modes  were  matched  within  2%  using  the  lengths  h = 
6.0  ft,  I2  = 5.5  ft  and  I3  = 6.0  ft.  From  this  process  the  model  torsional  stiffness’  kj2  = 3.6(106)  lb 
ft/rad  and  k23  = 1.69(106)  lb  ft/rad.  were  obtained 

Dynamic  analyses,  [2],  [3]  have  successfully  modeled  tube  response  using  proportional 
damping,  i.e.  [c]  = P[k].  In  the  case  of  the  XM291  P = 0.0015  sec  has  been  found  to  be 
reasonable.  A reasonable  estimate  for  trunnion  damping  is  clp  = 750  lb  ft  s/rad. 

The  first  attempt  to  determine  the  feedforward  transfer  function  using  Equations  (12)  and 
(13)  failed  because  the  some  of  the  roots  of  B4!  were  positive.  This  difficulty  was  eliminated  by 
using  the  alternative  strategy  described  above,  with  the  following  results. 

r44  = 1 + S2,  G(.  G22  = 0 (18) 

G<Gu=^r  (19a) 

G<Gn-^r  (19b) 

**21 

Gr  G\3  - ~~  (19c) 

**2l 
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The  det  [a]  plays  no  role  in  these  functions  because  it  is  canceled  by  ratioing  the  B’s.  Figures  3 
and  4 show  the  Bode  plots  of  GcGn  and  GcGi2-  The  transfer  function  GcGi3  is  zero  so  that  0P  is 
not  fed  forward.  The  numerators  and  denominators  are  all  fifth  order  polynomials,  so  that  the 
high  and  low  frequency  gains  are  bounded. 

The  feedback  transfer  function,  GCG22,  is  shown  in  Figure  5.  The  remaining  feedback 
transfer  function,  GCG24  plays  no  role  in  disturbance  rejection  in  this  case. 

The  elements  of  d22,  d23,  d32,  and  d33  of  [d]  were  found  to  be  identically  zero.  The 
remaining  no-zero  elements  of  Equation  (lOd)  and  (16)  yield  the  following  response  equations: 
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Figures  6 through  10  show  the  transfer  functions  required  by  the  equations  above. 

Figures  9 and  10  show  that  affect  of  the  trunnion  acceleration  on  0i  and  02  is  highly  attenuated  so 
that  large  angular  displacements  of  the  tube  are  not  required  to  achieve  stabilization. 

Figures  6 through  8 are  quite  similar.  It  appears  that  the  required  x c will  depend  largely 

on  the  trunnion  acceleration  and  pitch  rate  at  very  low  frequencies.  There  is  a considerable 
attenuation  of  the  disturbance  inputs  up  to  100  rad/sec  (~15  Hz)  with  a return  to  the  low 
frequency  levels  at  103  rad/sec  (160  Hz). 

CONCLUSIONS 

It  appears  that  the  results  previously  obtained  with  the  two-segment  model,  [1],  led  to  the 
erroneous  conclusion  that  the  effects  of  the  disturbance  could  not  be  entirely  repeated  from  the 
muzzle  angular  displacement.  The  analysis  of  the  three-segment  model  presented  suggests  that 
this  is  possible,  at  least  theoretically.  Of  course  the  unanswered  question  is  "what  would  be 
revealed  by  a higher  order  multi-segmented  model,  and  how  many  segments  are  enough." 

On  the  practical  side,  it  is  certain  that  the  transfer  functions  cannot  be  duplicated 
precisely.  There  are  four  of  these  that  must  be  implemented  with  reasonable  fidelity  to  achieve 
the  predicted  results  of  the  three-segment  model.  That  number,  along  with  their  input  signals, 
indicates  the  magnitude  of  the  task.  While  feedforward  and  feedback  control  has  long  been  used 
in  fire  control  it  is  hoped  that  this  paper  provides  some  guidance  in  their  use  when  tube  flexure  is 
a consideration. 

NOMENCLATURE 

fa]  - dynamic  matrix 
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[B]  = [C]  [I] 


[c]  = damping  matrix 

Ci 2,  C23  = damping  coefficients 
c ip  = trunnion  viscous  friction  coefficient 
[C]  = cofactor  matrix  of  [a] 

[d]  = disturbance  input  matrix 
{D}  = disturbance  vector 

Fc  = elevation  actuating  force 
Gc  = scalar  gain 

Gcd  = feedforward  transfer  function  vector 
Gd  = disturbance  transfer  function 
Gp  = plant  transfer  function 

G|  1,  G|i,  G 1 3 = feedforward  transfer  function  vector  components 
G22,  G24  = feedback  transfer  function  vector  components 
H = feedback  transfer  function  vector 

[l]  = forcing  function  matrix 
[k]  = stiffness  matrix 

k 1 2,  k23  = stiffness  coefficients 
kd  = drive  line  stiffness 
1 1, 12, 13  = segment  lengths 

[m]  = mass  matrix 

[q]  = disturbance  response  matrix 
xc  = elevation  actuator  displacement 
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X,p  = distance  from  trunnions  to  drive 
yt  = vertical  displacement  of  the  trunnion 
[r]  = cofactor  matrix  of  [q] 

P = proportional  damping  coefficient 
q i , r|2,  = center  of  mass  coordinates 
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m23 


*23 

c23 

*23 

C23  _ 

Aijj  + k22 

~k23 

~k23 

k23 

a 


11 


^21 

0 

0 


a,- 

0 

0 

12 

au 

-4c 

II 

Cl  22 

*23 

a2A 

«12 

— *21 

= ~kdX 

*32 

*33 

«34 

2 

a42 

a43 

«44_ 

fl24 

= «42 

= mX3S 

a22  m\\S  "*“(CI2  C\p)S  ■*"  i^d^tp  + kn) 


a23  ~ a22  ~ m\2S  C12S  k\2 
^33  ^ 22 ^ (^12  ^'23)*^  (^"12  ^“23) 

^34  ^43  ^23^  ^23*^  23 

«44  — W335  ^ C23‘5  ^23 
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[/]= 


•'ll 

0 

0 

0 


0 0/ 


1 22 


‘32 


‘42 


•“  23 
0 


14 


/ 


0 0 


/„=1.0 


1 22  ~ hm2 


1 23  ~ •'lp 

*24  = W 


/32  — (j]2m2  +^2^3) 

1,2  = ^3m3 


APPENDIX  D 


C„ 

C„ 

C,3 

[<=■]= 

C2] 

C22 

r 

'-'23 

Ql 

C32 

Q3 

C41 

c 

‘-'42 

c43 

C„  = 

* 22 

* 33*44 

+ la 

C22  = 

: *ti 

(*  33  * 44 

- a 

C33  = 

:*u 

(*  22  * 44 

- a 

= 

: *>1 

22  ^ 33 

- a 

C,2  = 

-C„ 

= “*2> 

(*33 

C13  = 

C„ 

= *2l(*33* 

CM  = 

Q, 

= ~*21 

(*  32  1 

C23  = 

C32 

= ~*u 

(*  32  1 

C24  = 

c 

^42 

- au(a32a43 

C34- 

r ; 

^43 

= -*u(* 

'22*43 

c2, 

c,. 


24  * 33  * 42  ^ 43  ^ 34  ^ 22  * 23  * 32  * 44  ) 


‘32  . 


det[*]  = *n  (*22*33*44  + *23*34*42  + ana„a24) 

~ *ll(*24*33*42  + *34*43*22  + *23*32*44  ) ~ *12*21  (*33*44  — *34*43  ) 


APPENDIX  E 


w- 


Bn 

Bu 

Bu 

Bu 

Bn 

B22 

Bn 

Bn 

Bn 

B22 

B„ 

Bu 

B,x 

B,2 

Bn 

B ,4 
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Bn  — Gu^n  / det [a ] 

Bn  =(C2,/22  +C31/32  +C4I/42)/det[a] 

Bn  =(C2i^23)/det^] 

BU  =(cn/i4  +C:i/24)/det[a] 


B2i  + (C12/n)/det[a] 

B22  = (C22/22  + C32/32  + C42/42  )/det[a] 


523  =(C22/23  )/det[fl] 

524  =(CnI»  +C22I2,)/det[a\ 


B3\  =(C13/u)/det[a] 

Bn  = (C23I22  + C33/32  + C43/42  )/det[a] 


533  =(C23723)/detH 


534  =(C13714  +C23/24  )/det['3] 


54!  =(Ci4/,i)/ciet[a] 

Bu  =(CiJi2  +C2J22  +C44/42  )/det[a] 


543  =(Q4A3)/det[a] 

544  =(C14/14+C24/24  )/det[a] 


APPENDIX  F 

Hr 


^ii  0 r41 

0 r22  0 f42 

0 '23  '*33  ^43 

0 r24  0 r44 


ru  = det[a]=  1 + B2iGcG22  + BAXGCG2 


G\  = -BxxGcG  22 

r22  = ^ + BAlGcG2A 

r22  = ~B2XGcG21 

>24  = ~B*lGcG22 


r4I  = ~B[\GcG24 
rn  = ~B2lGc  G24 
r43  = ~B3lGcG24 

^4  — 1 + B21GcG22 
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M= 


M2 


‘32 

^42 


i/n  t/12  af13 

^21  ^22  ^23 

dn  ~ B\2  BuGeGn 

^31  ^32  ^33 

d2x  = B22  — B2XGcGn 

_^41  ^42  ^43. 

dn  = 532  -B3lGeGu 

Bn-BuGcGr. 

dA\  = B42  -B4lGcGn 

B2i  - B2lGcGi2 

dn  = B\4  ~ BnGcGn 

B33  -B»GtGn 

dn  = B24  — B2XGcGu 

B4  3 — B4iGcGn 

dx  = 534  ~B3lGcGu 

d4i  = B44  — B4lGcGn 

APPENDIX  G 

zzt,  = cr.,  = /?rsf  fwo  natural  frequencies 

*23  ~ [(^j  + ®2  \m22m33  ~ ,fl2)m}2  ) 2^*23  + W22  )]^23 

+ W33  (t&i&l  )(w22ff233  - /M-j/Hj,  )/(/«j3  + 2m23  + W322  ) = 0 

*12  = (CT12  + ^22X/M22W33  - W:j#nj,)/»ljj  - (»!„  + 2w23  + mil)k2J  m» 
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